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ABSTRACT 

We present radio maps from the Very Large Array (VLA) for 16 sources detected in a sub-millimeter (submm) 
survey of the distant Universe. Our deep VLA 1 .4-GHz maps allow us to identify radio counterparts or place 
stringent limits (< 20/iJy in the source plane) on the radio flux of the submm sources. We compare the spectral 
indices of our sources between 850 /im and 1.4 GHz to empirical and theoretical models for distant starburst 
galaxies and active galactic nuclei (AGN) as a function of redshift. In this way we can derive redshift limits for 
the submm sources, even in the absence of an optical or near-infrared counterpart. We conclude that the submm 
population brighter than ~ 1 mJy has a median redshift of at least <z>^ 2, more probably <z>~ 2.5-3, with 
almost all galaxies at z 3> 1 ■ This estimate is a strong lower limit as both misidentification of the radio counterparts 
and non-thermal emission from an AGN will bias our redshift estimates to lower values. The high median redshift 
means that the submm population, if predominately powered by starbursts, contributes a substantial fraction of 
the total star formation density at high redshifts. A comparison of the spectral index limits with spectroscopic 
redshifts for proposed optical counterparts to individual submm galaxies suggests that half of the submm sources 
remain unidentified and thus their counterparts must be fainter than / > 24. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — infrared: galaxies 
— radio: galaxies 



1. INTRODUCTION 

Faint submm sources are likely to be highly obscured star- 
burst galaxies and AGN, within which optical/UV radiation 
from massive stars or an active nucleus is absorbed by dust 
and reradiated in the far-infrared. The dust emission peaks 
at A ~ 100/im and thus long-wavelength observations of dis- 
tant dusty galaxies can benefit as this peak is redshifted into 
their window. At 850 /im, this increase balances the geomet- 
rical dimming at higher redshifts, resulting in a constant flux 
density out to z ~ 5-10 and the opportunity to select very high 
redshift galaxies. 

A number of deep submm surveys have been published 
(Smail, Ivison & Blain 1997; Bargeret al. 1998, 1999b; Hughes 
et al. 1998; Blain et al. 1999a; Eales et al. 1999) providing 
counts which are in good agreement at 850-/im flux densities 
above 2mJy (the confusion limit of the blank-field surveys). 
The surface density of submm galaxies reaches ~ 2 per sq. ar- 
cmin by 1 mJy (Blain et al. 1999a). If these galaxies lie at z > 1, 
then they have bolometric luminosities of > 10 12 L Q and they are 
the distant analogs of the local ultraluminous infrared galaxy 
(ULIRG) population. However, the observed surface density of 
submm galaxies is several orders of magnitude greater than that 
expected from the local ULIRG population (Smail et al. 1997) 
indicating very substantial evolution of these systems in the dis- 
tant Universe. The integrated emission from this population can 
account for the bulk of the extragalactic background detected at 
850 /im by COBE (e.g. Fixsen et al. 1998), and hence confirms 
these galaxies as an important source of radiation in the Uni- 
verse (Blain et al. 1999a, 1999b). 
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To identify the era of obscured emission in the Universe, 
whether from AGN or starbursts, we have to measure the red- 
shifts of a complete sample of submm galaxies. Several groups 
have attempted this (Hughes et al. 1998; Barger et al. 1999a; 
Lilly et al. 1999). Hughes et al. (1998) concluded that the bulk 
of the population is at z ~ 2-4, based on photometric redshift 
limits for the probable counterparts of five submm sources in 
the Hubble Deep Field (HDF, c.f. Richards 1999 and Downes 
et al. 1999). Barger et al. (1999a) undertook a spectroscopic 
survey of the same submm sample analysed here and concluded 
that the median redshift was < z >^ 1 .5-2, with the bulk of the 
population having z ~ 1-3. Finally, Lilly et al. (1999) used 
archival spectroscopy and broad-band photometry of submm 
sources from the Eales et al. (1999) survey to claim that the 
population spans z = 0.1-3, with a third at z < 1. The differ- 
ences between these studies are significant and important for 
our understanding of the nature of submm galaxies. 

It is very difficult to achieve high completeness in optical 
spectroscopic surveys of submm galaxies (e.g. Barger et al. 
1999a) due to the very different behaviour of the K correc- 
tions for distant galaxies between submm and optical pass- 
bands. However, even a crude estimate of the median redshift 
of a complete sample of submm galaxies would provide a pow- 
erful insight into the relative dominance of obscured and unob- 
scured emission at different epochs (Blain et al. 1999b). 

In a recent paper, Carilli & Yun (1999, CY) demonstrated 
that using the spectral index between the submm (850 pm) and 
radio (1.4 GHz) wavebands, a^ 5 ® = 0.421og 10 (Sg5o/Si.4), it was 
possible to obtain crude redshift limits for distant dusty galax- 
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Radio Observations of Submm Galaxies 




Fig. 1. 15 X 15 aicsec VLA maps of the 16 submm sources for which we have radio observations. These are ordered as in Table 1 and the maps span a range in sensitivity and resolutions 
(§2.1). The panels are centered on the nominal submm position in each case and we show the typical error boxes for the sources. We mark the radio counterpart by a + where they are 
identified. The panels have north top and east left and correspond to lOO/i^ kpc at z > 1. The contour levels are at apparent 1.4-GHz flux densities of 3, 5, 10, 20, 50, lOOx the map noise 
listed in §2.1 (except for SMM J21536+1741 and SMMI14010+0252 where they are lOOx these values). 



ies, irrespective of the nature of the emission mechanism, AGN 
or starburst. CY employed a number of theoretical and em- 
pirical spectral energy distributions (SEDs) to investigate the 
range in a* 5 ° for different assumed SEDs and showed that these 
models adequately described the small sample of high-redshift 
galaxies for which both radio and submm observations were 
available. As pointed out by Blain (1999, B99), if we adopt 
lower dust temperatures for the submm population than are seen 
in the local sources used in CY's models, then the allowed range 
of redshifts is slightly lower for a given value of . Never- 
theless, the modest scatter between the models in CY suggests 
that this technique can provide useful limits on the redshifts of 
submm galaxies in the absence of an optical counterpart. 

In this paper we apply the CY analysis to deep radio obser- 
vations of a complete sample of submm galaxies selected from 
the SCUBA Cluster Lens Survey (Smail et al. 1998). Our aim 
is to constrain the redshift distribution of this population and 
in the process test the optical identifications and spectroscopic 
redshifts from Smail et al. (1998) and Barger et al. (1999a). We 
present the observations and their analysis in §2, discuss our 
results in §3 and give our main conclusions in §4. 

2. OBSERVATIONS, REDUCTION AND ANALYSIS 

The 850-/xm maps on which our survey is based were ob- 
tained using the long-wavelength array of the Sub-millimeter 
Common-User Bolometer Array (SCUBA, Holland et al. 1999) 
on the James Clerk Maxwell Telescope (JCMT) 8 . The details of 
the observations, their reduction and analysis are given in Smail 
et al. (1997, 1998) and Ivison et al. (1998). Each field covers 
an area of 5.2 arcmin 2 with a typical la sensitivity of ~1.7 mJy, 
giving a total survey area of 0.01 deg 2 . The median amplifica- 
tion by the cluster lenses for background sources detected in 
our fields is expected to be 2.5*1 5 (Blain et al. 1999a; Barger 
et al. 1999a), and so we have effectively surveyed an area of 
about 15 arcmin 2 in the source plane to an equivalent la sensi- 
tivity of 0.7 mJy. The follow-up of these submm sources also 
benefits from the achromatic amplification, which boosts the 
apparent brightness of counterparts in all other wavebands. 

2.1. Radio Observations 

All the radio maps used in this work were obtained with 
the VLA 9 at 1.4 GHz in A or B configuration, giving effec- 
tive resolutions of 1.5" and 5", respectively. More details of 
the reduction and analysis of these maps are given in the fol- 



lowing references (we list the VLA configuration and la map 
noise for each cluster): Morrison et al. (1999) for CI 0024+16 
(B/15 /iJy), A 370 (B/10 ^Jy) and CI 0939+47 (B/9 ,uJy); Ivison 
et al. (1999) for A 1835 (B/16 ^Jy) and Ivison et al. (2000) for 
MS 0440+02 (A/15 /zJy) and CI 2244-02 (A/17 /xJy). No deep 
radio map is available for A 2390, although shallower observa- 
tions were used to study the submm/radio spectral index of the 
central cluster galaxy (Edge et al. 1999). 

Radio counterparts were searched for around the nominal po- 
sitions of the submm sources based on the SCUBA astrome- 
try and 2a error-boxes of 6" or 8" diameter depending upon 
whether the submm source was respectively a 4a or 3a de- 
tection (Fig. 1). The size of these error-boxes includes both 
the systematic and random errors in the source positions and 
they have been confirmed as realistic using two sources with 
CO interferometric observations (Frayer et al. 1998, 1999). 
To assign radio fluxes to the individual submm sources we 
have adopted the conservative approach of identifying the ra- 
dio counterpart as the brightest radio source within the submm 
error-box (Fig. 1). The behaviour of the radio-submm spec- 
tral index is such that a brighter radio counterpart leads to a 
lower redshift estimate. Thus by selecting the brightest avail- 
able radio source, we should obtain a lower bound on the red- 
shift of the submm source. Apparent radio fluxes or 3a limits 
are listed in Table 1. Note that a faint radio source lies just 
outside the error-box of SMM J00265+1710, and so we adopt 
the radio flux of that source as a lower limit in this case. For 
SMMJ02399-0134, the radio counterpart is confused due to 
emission from a nearby bright cluster elliptical (Fig. 1), and so 
we only quote an approximate 1 .4-GHz flux for this galaxy. 

2.2. Radio-Submm Spectral Indices 

The 850-/im and 1.4-GHz fluxes or limits for 16 of the 
sources in Smail et al. (1998), for which we have radio obser- 
vations, are listed in Table 1 in order of their apparent submm 
fluxes, along with their proposed spectroscopic redshifts, z spec , 
from Barger et al. (1999a). Where a radio counterpart is iden- 
tified the spectroscopic redshift of the closest optical candidate 
is listed in the table. The errors on a\ 5 ^ are calculated assum- 
ing the la flux uncertainties in each band (§2. 1 and Smail et al. 
1998). For non-detections at 1.4 GHz we use the 3a flux limit 
of the relevant radio map. The two central cluster galaxies in 
our sample are not included in our analysis. 

Using the af 5 ® values or limits, the redshift ranges, z a , are 
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Table 1 

Spectral Properties of the Submm Sample 



Submm 
Source 


^850 

(mjy) 


Si .4 






Reliable 


4pec 

Uncertain 


Comments 


4cr Detections 
















SMMJ02399-0136 


25.4 


526 


0.70 ±0.02 


1.1-2.9 


2.80 




L1/L2 - Seyfert 2 merger" 


SMMJ00266+1708 


18.6 


100 


0.95 ±0.04 


2.0-4.5 




0.44 


M8 - l" away from radio source 


SMM J09429+4658 


17.2 


32 


1.14 ±0.07 


>3.9 






H5 - ERO* 


SMM J 14009+0252 


14.5 


529 


0.60 ±0.03 


0.7-2.3 






J5 C - no spectroscopic observation 


SMM J 1401 1+0252 


12.3 


115 


0.85 ±0.05 


1.7-3.8 


2.55 




J1/J2 - starburst merger^ 


SMM J02399-0134 


11.0 


~500 


0.55 ±0.04 


0.6-2.1 


1.06 




L3 - Seyfert 1.5/2 ring galaxy* 1 


SMM J2247 1-0206 


9.2 


<65 


^ n on 4- 08 


>1.8 




1.16 


P4 - weak AGN? 


SMM J02400-0134 


7.6 


<33 


> 0.99±0.11 


>2.4 






Optical blank field 7 


SMM J04431+0210 


7.2 


<70 


> 0.84 ±0.10 


>1.6 






N4 - ERO* 


3cr Detections 
















SMMJ21536+1742 


6.7 










1.60? 


No deep radio map 


SMM J00265+1710 


6.1 


<110 


> 0.73 ±0.07 


>1.2 




0.21 


M6/M10 - radio ID outside error-box 


SMM J22472-0206 


6.1 


<50 


> 0.87 ±0.12 


>1.7 




2.11? 


P2 


SMM J00266+1710 


5.9 


<33 


> 0.93 ±0.12 


>2.9 




0.94 


M3 - faint radio source at map limit? 


SMMJ00267+1709 


5.0 


<30 


> 0.92±0.13 


>1.9 






Optical blank field 7 


SMM J04433+0210 


4.5 


70 


0.75 ±0.12 


1.3-3.2 






No spectroscopic observation 


Cluster Galaxies 
















SMMJ21536+1741 


9.1 


226 x 10 3 


-0.58 ±0.05 




0.23 




cD - central galaxy in A 2390* 


SMM 114010+0252 


5.4 


28.6 x 10 3 


-0.30 ±0.06 




0.25 




cD - central galaxy in A 1835 s 



REFERENCES.— a) Ivison et al. (1998) - b) Smail et al. (1999) - c) Ivison et al. (1999) - d) Barger et al. (1999a) - e) Soucail et al. (1999) -/) Smail et al. (1998) — g) 
Edgeetal. (1999). 



derived from the extremes of the predictions from the four CY 
models (two empirical SEDs representing Arp 220 and M 82, 
and two models with dust temperatures of ~ 50-60 K and 
emissivities of (3 = 1.0 or f3 = 1.5) and a further model from 
B99, with f3 = 1 .5 and a = 30 K to illustrate the minimum 
possible redshift assuming a very low r<j. For a model SED 
at z > 1, a variation of 6(3 = +0.2 or 5T^ = -10K results in an 
change of 5a\ 5 ^ < +0.1, equivalent to an uncertainty in the de- 
rived redshift of Sz/z ~ 0.2. 

There are three caveats to bear in mind when using a\ 5 ® to 
estimate redshifts for distant galaxies. First, most of the dis- 
tant galaxies which CY used to compare with their model pre- 
dictions show some signs of AGN activity. If these AGN also 
contribute to the 1 .4-GHz non-thermal emission of the galaxy 
they will lower the observed aj 5 4 ° values (as some obviously do 
in Fig. la of CY). This will mean that any radio-quiet submm 
sources could lie at the high end of the predicted a* 5 ® range 
at each epoch. Secondly, the CY and B99 models we use as- 
sume effective dust temperatures for the galaxies, > 30 K. If 
the dust in distant obscured galaxies is much cooler then this, it 
will again shift the predicted redshifts systematically lower (see 
B99). Finally, as mentioned in CY, the effects of inverse Comp- 
ton scattering of radio photons off the microwave background 
may reduce the radio luminosities of star-forming galaxies at 
z > 3 and hence increase a\ 5 ® for the most distant galaxies. 
Nevertheless, the relatively good agreement shown in CY of 
the spectral indices of distant galaxies with the models is an 
important confirmation that a\ 5 ® indices can be used to derive 
robust lower limits to the redshifts of submm sources without 
reliable spectroscopic identifications. 

3. RESULTS AND DISCUSSION 

We show in Fig. 2 three cumulative redshift distributions for 
the population representing extreme interpretations of the z a 
limits from the spectral index models. We see that even mak- 
ing the most conservative assumptions about the likely redshifts 



from the cq f indices we still predict a median redshift for the 
submm population above an intrinsic 850-/jm flux of 1 mjy of 
<z> ~ 2, and more likely closer to z ~ 2.5-3. 

Comparing the cumulative redshift distribution to that de- 
rived from the (incomplete) spectroscopic study of this sample 
(Barger et al. 1999a) we see broad similarities. However, com- 
parison of redshifts for individual sources from the two stud- 
ies (Table 1), while showing good agreement for those submm 
sources with reliable identifications (e.g. Ivison et al. 1998, 
1999) also indicates that the majority of the uncertain spec- 
troscopic IDs are likely to be incorrect. Barger et al. (1999a) 
obtained spectroscopy of most of the possible optical counter- 
parts within each submm error-box. We can therefore state that 
the true submm sources must be fainter than the faintest spec- 
troscopic target. Including the two optical blank-field sources 
already known (Smail et al. 1998), we conclude that approx- 
imately half of the submm population are therefore currently 
unidentified. These submm sources have no radio counterparts 
and are too faint for optical spectroscopy, / > 24, their identifi- 
cation will thus be very difficult. 

Our median redshift is compatible with the results of Hughes 
et al. (1998) and CY for the five submm sources in the HDF 
based on analyses of their SEDs and radio-submm indices. The 
only other submm survey for which spectroscopic redshift in- 
formation has been published is by Lilly et al. (1999) for the 
Eales et al. (1999) sample — they find <z> ~2 — similar to 
our median redshift. However, Lilly et al. (1999) claim that a 
third of the submm population lies at z < 1 ; in contrast, we find 
no galaxies in our field sample at z < 1 . This apparent contrad- 
ication may result simply from the small sizes of the current 
samples or might indicate that foreground bright optical galax- 
ies are lensing the distant submm sources detected in the field 
surveys (see Blain, Moller & Mailer 1999; Hughes et al. 1998). 

The detection rate of radio counterparts to the submm 
sources is higher for the intrinsically brighter sources. All the 
submm sources with observed fluxes above ~ lOmJy (intrinsic 
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fluxes of > 4mJy ) have radio counterparts, while the majority 
of the fainter sources do not (this is consistent with CY's re- 
sults in the HDF). The detections and astro me try of the fainter 
sources are sufficiently reliable that this result is not due to spu- 
rious detections (see Ivison et al. 1999). Making the conser- 
vative assumption of placing all of the non-detections at their 
lower bounds on a\ 5 ®, we find that this distribution is consis- 
tent with being drawn from the aj 5 " distribution of the brighter 
sources with a probability of P ~ 0.4. However, simply com- 
paring the aj 5 ^ indices for the two subsamples we see that half 
of the radio-detected bright submm sources have aj 5 4 ° values 
lower than the lowest limit on the undetected sources, the like- 
lihood that this occurs by chance is only f~3x 10~ 4 suggest- 
ing that there may be real differences between the aj 5 4 ° val- 
ues for the two subsamples. Several factors could cause this, 
most simply the apparently fainter submm sources may be at 
higher redshifts (a correlation which could naturally exist in a 
low density Universe). Alternatively, intrinsic differences in the 
spectral indices of fainter sources would occur if they contain a 
lower fraction of radio-loud AGN or have typically cooler dust 
temperatures (or higher emissivities). Both effects are plausi- 
ble given what we know about the correlations of AGN frac- 
tions, dust temperature and emissivity with luminosity in local 
ULIRG samples (Sanders & Mirabel 1996). Further detailed 
observations of both distant and local ULIRGs are needed to 
distinguish between these possibilities. 




CY+B99 minimum-z 

CY only minimum-z 

CY+B99 flat P(z) 

Borger et al. (1999) 
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Fig. 2. The cumulative redshift distribution for the full submm sample. We have used 
the spectroscopic redshifts of those sources thought to be reliable (Table 1) and combined 
these with the probable redshift ranges of the remaining sources derived from their aj 5 4 ° 
indices or limits. The solid line shows the cumulative distribution if we assume the mini- 
mum redshift distribution which is obtained if all sources are assumed to lie at their lower 



Z a limit given in Table 1 (the dashed line is the equivalent analysis but restricted to just 
the CY models). The effect of non-thermal radio emission, which drives down the a] 5 ° in- 
dices, means that this is a very conservative assumption if some fraction of the population 
harbor radio-loud AGN. The dot-dashed line assumes a flat probability distribution for the 
sources within their z a ranges and a maximum redshift of z = 6 for those sources where 
we only have a lower limit on a^ 5 ® . Finally, the dotted line is the cumulative redshift dis- 
tribution from Barger et al. (1999a) with two of the source identifications corrected as in 
Smail et al. (1999) and all blank-field/ERO candidates placed at z = 4. 

The relatively high median redshift we find for the submm 
population, <z> <~ 2-3, indicates that their equivalent star for- 
mation density at these epochs is around 0.5 M Q yr" 1 Mpc" 3 
(Blain et al. 1999a), roughly three times that seen in UV- 
selected samples (Steidel et al. 1999). Emission from dust 
heated by obscured AGN will reduce this estimate, but it is diffi- 
cult not to conclude that the submm galaxies contain a substan- 
tial fraction of the star formation in the high redshift Universe. 

4. CONCLUSIONS 

• We present radio maps of 16 galaxies selected in a deep 
submm survey. We combine submm and radio fluxes (or limits) 
to determine the radio-submm spectral indices of these galaxies 
and interpret these using model predictions to derive the red- 
shifts for a complete sample of faint submm galaxies. 

• We find a median redshift < z > ~ 2 for the submm popu- 
lation down to #850 ~ 1 m Jy under conservative assumptions, 
and <z> <~ 2.5-3 for more reasonable assumptions. Median 
redshifts below < z > -C 2 are only possible if the bulk of the 
emission is coming from dust at <C 30 K (compared to the 
40-50 K typically seen in well-studied, distant submm galax- 
ies, or their low-redshift analogs: ULIRGs). As a result we 
find no evidence for a significant low-redshift, z < 1 , tail in our 
distribution in contrast to Lilly et al. (1999). 

• We compare the individual redshifts estimated from a \ 5 ® with 
the spectroscopic observations of proposed optical counterparts 
of the submm sources. We find that the majority of the 'uncer- 
tain' spectroscopic identifications from Barger et al. (1999a) are 
likely to be incorrect. We conclude that the true counterparts lie 
at higher redshifts and are intrinsically very faint, / > 24, mak- 
ing the prospects for a complete optical spectroscopic survey of 
the submm population bleak. 
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